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The equilibrium size and temperature limits of thermally and chemically equilibrated e + e~7 
plasma drops are investigated at a given energy content. For a plasma to be equilibrated it must be 
opaque to electron and photon interactions. The opaqueness condition is determined by comparing 
plasma size with the mean free electron and photon paths. We calculate those paths using thermal 
Lorentz-invariant reaction rates for pair production and electron (positron) and photon scattering. 
The range of the corresponding plasma temperature and size is evaluated numerically. Considering 
the energy and size we find that the opaque and equilibrated plasma drop may be experimentally 
attainable. 
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I. INTRODUCTION 

In [l[ we considered relativistic EP 3 (e~, e + , 7 plasma, 
electron, positron photon plasma) in equilibrium and 
have studied heavy (muon, pion) particle production 
rates. Here we extend considerably these methods and 
study the production and equilibration of the compo- 
nents e~,e + ,7 of the EP 3 plasma in order to explore 
the constraints on how such a plasma can be formed, 
and when it falls apart into free-streaming constituents 
(freeze-out condition). We evaluate constraints arising 
from the requirement that a plasma drop ought to be 
opaque to constituent particle scattering, and include in 
this consideration the constraint on size and tempera- 
ture arising from laser pulse energy available to form the 
plasma. Furthermore, we can determine under what con- 
ditions the plasma drop will not only reach thermal (spec- 
tral) but also chemical (yield of particles) equilibrium. 

In laboratory experiments involving, e.g. ultra intense 
laser pulses to form strong fields, electron-positron pairs 
are produced; for review see 0t1|- They rescatter or an- 
nihilate and produce hard photons. In a large enough do- 
main the e + e~ pairs and hard photons 7 will equilibrate 
thermally and chemically. This kinetic equilibration pro- 
cess cannot be resolved in the general study we present 
here. Instead, we focus our attention on the properties 
of the e + e~7 plasma drop near equilibrium in its final 
evolution stage. We are assuming EP 3 was formed in the 
vacuum state (particle-antiparticle symmetry) and keep- 
ing the total zero spin. 

In laboratory experiments an initially dense plasma 
fireball with radius R p \ will dilute, expanding at plasma 
speed of "sound" , nearly speed of light. The lifespan in 
a fast expanding plasma drop is 



—Rph 



We consider the plasma drop of finite size. However we 
neglect the energy loss from the radiation of photons from 
the plasma surface. We assume that if a plasma is opaque 



then the energy loss from a surface is insignificant during 
the plasma lifespan, which is defined by the laser pulse 
duration. If plasma is nonopaque, then energy loss oc- 
curs from the whole plasma volume, signaling the end of 
equilibrium and onset of transport regime not addressed 
here. 

There are three different conditions which we can track 
in the plasma state which are usually passed in sequence 
in time: 

1. At highest density (i.e. temperature T) we reach 
particle yield chemical equilibrium. This is the con- 
dition when the e + e~-pair annihilation into and 
production by photons occur at equal rate. 

2. At an intermediate density many scattering pro- 
cesses can occur, but not as many as needed to 
maintain chemical equilibrate. That is the opaque- 
ness domain. 

At this time also most photons will leave the plasma 
free-streaming 

3. At yet a lower density, freeze-out of e + ,e~ occurs; 
that is, these particles will leave the plasma without 
further rescattering. 

Let us look at how these conditions arc obtained, in re- 
verse order: for the particle freeze-out (free-streaming) 
condition we consider 



Li — Rr) 



free — streaming. 



(1) 

" , 7 which 

is computed below, and compare it to the plasma radius 
R p i . For plasma to be opaque we require 



Here we use the thermal free path Li 



Li < 



R p i 



opaque. 



(2) 



Here and below the factor 3 is a putative choice, made 
in view of other studies of kinetic plasma dynamics. The 
e + e~ pair and photon chemical equilibration condition 
is approximately the same as the opaqueness condition, 
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a) Pair production: 

e + 



Vs 





b) Compton scattering: 






FIG. 1: Feynman diagrams for (a) the creation of the electron- 
positron pair on top and (b) Compton scattering on the bot- 
tom, where incoming particles are on the left side. 



a) M0ller scattering: 




b) Bhaba scattering: 





FIG. 2: Feynman diagrams for (a) M0ller scattering on top 
and (b) Bhabha scattering on bottom, where incoming parti- 
cles are on the left side. 



Eq.([2]) for photons participating in the pair production 
reaction Eq.(|4]): 



L 77 _j. ee < chemical equilibrium, 

o 



(3) 



The free-streaming and opaqueness conditions are up to 
numerical evaluation the same. 

We consider the mean free path L 7 and 7 equilibration 
by Dirac pair production/annihilation reaction: 



7 + 7 O e + e 



and by Compton scattering: 



7 + e ± -H- 7 + e ± , 



(4) 



(5) 



for which the Feynman diagrams are shown in Fig. [T] 

However, the evaluation of chemical equilibration of 
electrons and positrons is in essence different, these par- 
ticle will continue to rescatter even after all photons have 
left, and e~, e + can pair annihilate into photons. We ob- 
tain the mean free path L e - e + and the e + e~ pair equi- 
libration considering M0ller: 



e ± e ± «-» e ± e ± , 



and Bhabha scattering: 



(6) 



(7) 



for which the Feynman diagrams are shown in Fig. [5] 

We evaluate using a Lorentz covariant procedure the 
invariant reaction rates using methods developed first 



for strange quark pair production in the quark gluon 
plasma 0, @. We extend these methods to scattering 
reactions which do not change the number of particles. 

Our kinetic equations describing the rates and parti- 
cle multiplicity evolution are obtained from the Boltz- 
mann equation, integrated over the momentum part of 
the phase space, assuming thermal distribution. As a 
result, we obtain thermally averaged reaction relaxation 
time. This approach is appropriate to identify the do- 
main of T in which reactions are fast enough to main- 
tain thermal and chemical equilibrium. However, a fu- 
ture study will need to address how the thermal distri- 
bution is established. We also see a need to consider in a 
more complete transport approach the motion of parti- 
cles found in high energy tail of distribution, which will 
be first to escape the plasma drop. However, these par- 
ticles contribute little to the bulk properties. 

We observe that the rates for reactions considered here 
had been obtained by the other methods, using cross sec- 
tions Q and applied for astrophysical environments 
Our results differ from those presented in [8| , which found 
that under certain conditions plasma equilibration occurs 
primarily by pair production, and the Compton scatter- 
ing has small effect. We show that the Compton scat- 
tering is always the dominant process near equilibrium. 
The work to reconcile [8| with our covariant results con- 
tinues. One possible source of discrepancy is that in Q 
the Lorentz-invariant form of distribution function [see 
Eq.®] is not used in the rate evaluation procedure which 
involves a later change of reference frames. Therefore our 
Lorentz covariant results are different and amend this 
earlier work. 
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We do not consider the three-body reactions. Their 
rates are proportional to an additional factor a/71", and 
thus in general are slower than two-body reactions con- 
sidered here. We find that Bhabha and M0ller processes 
are much faster than other reactions considered. How- 
ever, if Bremsstrahlung is about a = 1/137 times slower 
than Bhabha and M0ller scatterings, it could still be 
faster, or of the same order of magnitude, as Compton 
and pair production. Bremsstrahlung and the related 
backward processes could therefore somewhat shorten 
photon equilibration, perhaps visible in some domain of 
T. We leave this to future study. 

Our paper is organized as follows: In section [TT] we 
review the physical properties (subsection III Ap of EP 3 
plasma, and present (subsection III B|) the master equa- 
tions for chemical equilibration of the plasma components 
under the assumption that particles are in thermal equi- 
librium. We introduce r™, the relaxation time for photon 
chemical equilibration, in terms of thermal reaction rates, 
and relate crrt to the mean free path for pair production. 

In Sec. IIIII we obtain the analytical expressions for the 
thermal mean free paths for scattering in EP 3 plasma. 
In Sec. 1111 Al we review equations for mean free paths 
and Lorentz-invariant rate calculations. In Sec. 1111 Bl 
we review cross sections for pair production and anni- 
hilation and Compton scattering as functions of the en- 
ergy of interacting particles. In Sec. IIII CI we obtain the 
Lorentz-invariant rate for pair production and annihila- 
tion. Similar invariant rate calculations were done before 
for the strange quark pair production in Refs. @ and @. 
We perform the required integration by a method simi- 
lar to Q, and muon pair production [l[ we considered 
earlier. In the following sections we extend that method 
to Compton (Sec. IHTD]) , M0ller and Bhabha (Sec. \UTE\i 
invariant rates. To remove the Coulomb divergence in in- 
tegrals we introduce a thermal screening photon (Debye) 
mass, induced by plasma oscillations. 

Sec. IIVI contains our numerical results. In Sec. IIV Al 
we look in depth at photon-electron interactions and the 
Thompson limit. In Sec. IIV Bl we evaluate the photon 
and electron scattering lengths, and discuss maximum 
possible plasma size at a given energy. This evaluation is 
carried out taking account of all quantum physics Fermi 
and Bose statistic effects in order to be applicable in the 
high density domain, and we apply the Lorentz-invariant 
rates method which include medium quantum effects. We 
connect our results to the Boltzmann limit in Sec. IIV CI 

In Sec. [V] we discuss the physical consequences, appli- 
cations, and present our conclusions. 



II. STATISTICAL PROPERTIES OF EP d 
PLASMA 



ativistic covariant Fermi and Bose momentum distribu- 
tions: 



fe± 



fl 



1 



(8) 

Tj (i = e, 7) is the fugacity of a given particle. When 
T» = 1, the particle yield is in chemical equilibrium. For 
T s ; = 1 this distribution maximizes the entropy content 
at a fixed particle energy Q. 

The Lorentz-invariant exponents involve the scalar 
product of the particle 4-momentum p^ with the local 
4- vector of velocity . In the absence of matter flow we 
have in the rest frame of the thermal bath (laboratory 
frame) 



(9) 



When the electron chemical potential v e is small, v e <C T, 
the number of particles and antiparticles is the same, 
n e - — n e + . Physically, it means that the number of 
e + e~ pairs produced dominates the residual matter elec- 
tron yield. Here we will set v e = , and will consider 
elsewhere the case for very low density plasma where 
chemical potential may become important. 

The particle density and plasma energy density can be 
evaluated using relativistic expressions: 

m = [ gifi(p)d 3 P ; E= [j29iEifi(p)d 3 pV, (10) 



where Ei — y/ mf + p 2 , fi (p) is the momentum distribu- 
tion of the particle i £ 7, e 1 , and <?, is its degeneracy; 
gi = 2 for photons, electrons and positrons. In Fig. [3] we 
present electron-positron and photon densities as func- 
tions of temperature. At our range of temperatures their 
values are much larger than atomic density, 10 23 cm~ 3 . 
After integration we obtain 



- = e = —g{T)T\ 
V 30^ ' 



(11) 



At a temperature T <C m e we only have truly mass- 
less photons and g(T) = 2 7 . Once the temperature ap- 
proaches and increases beyond m e , we find g ~ g'(T) ~ 
2 7 + (7/8)(2 e - +2 e+ ) = 5.5. 

In a classical case, T << 1, there is no difference be- 
tween Bose and Fermi particles. We have for massless 
particles (m/T ->• 0) 



E = 3NT, 



(12) 



where 5 = 6. 



Relativistic Gas 



B. Approach to equilibrium 



Up to QED interaction 0(a) effects we can describe 
the particle content in a plasma using, respectively, rel- 



In our approach new photons are produced by a pair 
annihilation reaction, Eq Q . If we assume that thermal 
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FIG. 3: Electron-positron (blue solid line) and photon (green 
dash-dotted line) densities as functions of temperature. 



FIG. 4: The relaxation time for pair and photon chemical 
equilibration. 



equilibrium is established faster than chemical, then the 
photon density evolution equations and chemical equilib- 
rium conditions arc similar to those for muon production, 
studied in 



1 diV 7 
V~dT 



= (Tg - Tl)R 



where 



R 



1 dW, 



1 dW P 



-f 2 

7 



dVdt 



T 2 



dVdt 



(14) 



dW 77 _ >e + e - jdVdt and dW e + e -^ rll /dVdt are Lorentz- 
invariant rates for pair production and annihilation re- 
actions, respectively. Rate i? 77< _> e + e - does not depend 



on T 

33 



in limit of a classical Boltzmann distribution. The 

< 1 



EP plasma is in relative chemical equilibrium for T 7 
when 



(15) 



In relative equilibrium at a given instant the considered 
reaction, for example (j4j, is in equilibrium. However 
for other reactions this condition may not be an equilib- 
rium (for example for reaction ee 4-> eej) and then the 
system can move to another relative equilibrium condi- 
tion with time. The relative chemical equilibrium oc- 
curs when crossing different time scale governs relevant 
"chemical" reactions. At T fii7 — > 1 (and for all others 
particles in the system, if any, Tj — > 1) all reactions 
are in equilibrium and the system achieves full chemical 
equilibrium. 



We introduce pair production relaxation time defined 



by: 



1 drij/dT 



■ 7 



2T e R 



■77-f^-e+e 



(16) 



(13) then for the simplest case T(t) = Const, and R(t) 



Const, the equation for T 7 is 



c£T 7 
~~dt 



(T 2 e - T*) 



2T P 



(17) 



The relaxation time vs temperature for photon and 
e + e~ pair chemical equilibration is shown in figure [4] 
The rates and relaxation times are discussed in depth 
in Q. 

To understand better the choice of definition of relax- 
ation time t 77 , Eq . (fT6|) . we assume T e = const, and in- 
troduce variable 7 = T 7 /T e , which shows deviation from 
chemical equilibrium. Then the equation for 7 is 



7 2 ) 



1 



ZT 77 



(18) 



The relaxation time is defined as the time for particle 
multiplicity to reach the magnitude of equilibrium value. 
Note that the relaxation time for a particle (in this ex- 
ample photon) production in a two-to-two-particles reac- 
tion, Eq. dTBl . changes by a factor T" 1 , where i denotes 
the initial particle in a reaction (here e^). The physical 
reason why the relaxation time is proportional to T^ 1 is 
that the collision rate drops by this factor, due to reduced 
density of plasma. 
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In the simple case considered here, e + e~ pair produc- 
tion and annihilation are chemically equilibrated when 
condition Eq. (|3j) is satisfied. 



III. MEAN FREE PATHS OF PHOTON AND 
ELECTRON (POSITRON) IN e+e - 7 PLASMA 



where I = 1, if there are identical particles in the initial 
or final state. 

Rate -R12.H.34 is the same for reactions in both direc- 
tions, because of time reversal symmetry of the matrix 
element: 



A. Reaction rate and mean free paths 

In order to be in thermal and chemical equilibrium 
plasma must be opaque for the reactions, which estab- 
lish this equilibrium. The major reactions which may 
establish thermal and/or chemical equilibrium between 
photons and e + e~ pairs are Compton scattering and pair 
production and annihilation. e + e~ can also participate 
in M0ller scattering Eq.© and Bhabha scattering Eq.Q. 

Here we consider the reaction 



1 + 2 o 3 + 4, 



(19) 



For reaction Eq. (fHj|) , the thermally averaged cross sec- 
tion is 



(fCTi2^34) 



Ill 2-0-12034 

nin 2 



(20) 



(velocity v = c for photons scattering, we take c = 1) 
and the mean free path is 



Li 



1 



(2) 



"1(2) 



"2(1) (WCT12^34) TiT 2 i?i2o 



34 



(21) 



where v is the relative velocity of interacting particles 
and (Ti2->.34 is cross section. 

Rate R is connected to Lorentz-invariant reaction rates 
by Eq.([l4]). The equation for rate R is [!, H] 



I? 



1 



12«34 — 



d 3 px f d 3 p 2 f d 3 p 3 f d 3 p 4 



1 + / (2tt) 8 J 1E X J 2E 2 J 2E 3 J 2E 4 

X(5 4 (pi + p 2 ~ P3 ~ Pi) ^ I (PlP2 |Mi2^34| PlPi) 1 2 X 



spin 



u- (P1+P2) 



\(piP2 \M 12 ^ 3 4\P3P4)\ 2 = \(P3P4 \M 34 ^ 12 \ Pl p 2 )\ 2 ■ (23) 



In the evaluation of matrix element we will use Mandel- 
stam variables: s, u, and t. For reaction (JTHJ) we have 



s = { Pl +p 2 ) 2 ; u=(p 3 -p 2 ) 2 ; t=(p 3 - Pl ) 2 ; (24) 



and 



t = mi 



If we exchange particles 2 and 3asl+3f>2 + 4 
or p 2 <-$■ p 3 , the matrix element for this reaction can be 
obtained from the matrix element for reaction Eq. (|19[) by 
crossing s and t: 



|M 12 ^ 34 (s, *, u)\ 2 = ± |Mi 3 ^2 4 (t, s, u)| 2 



(25) 



/i(pi)/2(p2)/3(P3)/4(P4)(T 1 T 2 T 3 T 4 )- i e— 



(22) 



The sign before the matrix element may change to op- 
posite depending if one of the particles, 3 or 2, changes 
from particle to antiparticle after crossing. 

For example, for pair production and annihilation and 
Compton scattering particles 1,2 are photons and parti- 
cles 3,4 are electron and positron. The matrix elements 
for these reactions is defined by Feynman diagrams in 
Figs. [U a) and (b), respectively. The matrix element av- 
eraged over final states spins for pair production and an- 
nihilation is [H, 



,, 99/ /to 2 m 2 \ 2 4m 2 4m 2 m 2 — u m 2 — t\ , , 

M 770ee 2 = 32ttV -4 -5— - + — + -5— + — + - + — , 26 

\ \m z — t m z — u J m z — t m z — u m z — t m z — u I 

The matrix element for Compton scattering is 

.,, ., „„ 9 9 / , / m 2 m 2 \ 2 4m 2 4m 2 m 2 — u m 2 — s\ . . 

M w 2 = 327rV 4 + — -— - - - , 27 

\ y m A — s — u J m A — s m z — u m z — a — u J 

The matrix elements for these reactions are connected by Eq. ([25]) . 
Similarly, crossing s and u allows us to find 

\M 12 ^ 34 {s,t,u)\ 2 =±\M 23 ^ 14 (u,t,s)\ 2 . (28) 
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The matrix elements for M0ller and Bhabha scattering (diagrams are in Figs. [5] (a) and (b), respectively) are connected 
by Eq. (|28|) . The M0ller scattering matrix element is (similar to [T3 |) 

2 2 j s 2 + u 2 + 8m 2 (t-m 2 ) s 2 + t 2 + 8m 2 (?i - m 2 ) (a - 2m 2 ) (s - 6m 2 ) ) 
1 e±e±l ~ \ 2(i-m 2 ) 2 + 2( U -m 2 ) 2 + (t - m 2 )(u - m 2 ) /' ^ 

To obtain the Bhabha scattering matrix element we need to cross u and s, Eq. (|2"5|) in M0ller scattering matrix element 
\. We obtain 

2 o6 2 2 | ^+ii 2 + 8m 2 (t-m 2 ) u 2 +f 2 + 8m 2 (.s-m 2 ) (u - 2m 2 ) (u - 6m 2 ) | 
|Me+e -' =2?ra \ 2(*-m 2 ) 2 + 2( S -m 2 ) 2 + (i - m 2 )(.s - m 2 ) /' (30) 

and the photon thermal (Debye) mass m 7 will be defined in Sec. IIIIEI 

B. Pair annihilation and production and Compton scattering cross sections 

In this section we consider cross sections for pair annihilation and production and Compton scattering as functions 
of energy. The cross section a(s) can be obtained by averaging the matrix element over the t variable. Similar 
calculations were done in [loj for heavy quarks production in quark gluon plasma. In the case of Compton scattering 
we obtain 

1 f° 

CT 7e ^ 7 e(s) = — — ^ / dt\M ±\ 2 

2na 2 { 2{s - m 2 ) 2 s - (s - m 2 f o „ / 8m 4 4 2 , 2 A, s\ 

' + 8m 2 - s + 4m 2 - (s - m 2 ) In — =■ (31) 



(s — m 2 ) 2 \ 2s 2 
For pair production and annihilation we have (similar results as in fill ]) 

1 



s(l-^/l-4m 2 /s)/2 
(l + ^/l-4m 2 /s)/2 



cr 77— >ee( s ) — -. „ 9 / ^|A^77V>ee| 



AlTOt 2 1(2 a 2 4\ + vs — 4m 2 



(s 2 + 4m 2 s - 8m 4 ) In v v / = - (s + 4m 2 ) y/ 3 (s - 4m 2 ) . (32) 
•y/s — v s — 4m 2 / 



7^m 2 -s(l-y/l-4m 2 /s)/2 
m 2 -s(l+^/l-4m 2 /s)/2 



2na 2 1/2 ,2 4\ + vs - 4m 



2 



s 2 + Am 2 s - 8m 4 ) In v = - (s + 4m 2 ) 0s (s - 4m 2 ) . (33) 



s 2 (s-m 2 ) \ yG-Vs- 4m 2 



Note that there is an extra 1/2 factor in the equation 
for pair annihilation because we have two identical par- 
ticles; that is a symmetrical wave function in the final 
state. Similarly, in the reaction rate R we add an addi- 
tional factor 1/2 when there are initial identical particles. 
In the backward reaction there is also a factor 1/2 from 
the definition of cross section. Therefore the rate is sym- 
metrical in both reaction directions: pair formation and 
annihilation. 



In Fig. [5]we show cross sections (|3"Tj) - (l3"2l as functions 
of total energy of particles in reaction s 1 / 2 (center of 
mass frame) for pair production (solid blue line) , annihi- 
lation (dash-dotted green line) and for Compton scatter- 
ing (dashed red line). 

At high energy Exot >> m the cross section depen- 
dence on energy is a power law, 

a = A(^-Y. (34) 
VlMeVy V ; 
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FIG. 5: The cross sections for pair production (solid blue line) 
and annihilation (dash-dotted green line) in the center of mass 
frame, and for Compton scattering (dashed red line) in the 
electron rest frame are shown as functions of total energy of 
interacting particles -Etot- 



TABLE I: Values of N and A for power law energy dependence 
of cross section high energy tail, Eq. (|34[) . for cross sections 
presented in Fig. \S\ 



Reaction 


N 


A (mb) 


7e «-> 7e 


- 1.7 


330 


77 o e + e~ 


-1.7 


410 


e + e~ -o 77 


-1.7 


309 



The coefficients A and N are shown in Table fl] 

The decrease of cross sections with particles energy 
increase could result in the early escape (freeze-out) of 
very high energy particles from the plasma which would 
impact heavier particle (muon, meson) production. How- 
ever, this has no material impact on our present consid- 
erations. 



q'=_p 4 +P35 p' = ^(P4~P3); 



(36) 



and then we have q = q$ — q = q' Q — q 12 = s > 2m 
and 



P3= 2 p; P4=2+P- 



(37) 



Using 



d 3 p 
2E 



d 4 p5{p 2 - m 2 )e{ Po ) = J d 4 P S Q (p 2 - m 2 ), 

(38) 

and p\ 4 — m — and p\ 2 = 0, we obtain: 



R 



77~f4-ee 



d 4 q / d 4 p / d 4 p'S (pI) So (pI) 



2(2tt) 8 

x<5 {pj - to 2 ) 6 (pj - m 2 ) T; 2 T- 2 exp (q /T) 

< Y, l^77-ee| 2 / 7 (p?) fe ( P ° 3 ) / 7 { P ° 2 ) fe (rf) • (39) 



The integrals from Eq. ([39|) can be evaluated in spher- 
ical coordinates. The angle coordinates are chosen with 
respect to the direction of ~q* =pt+pt: 

% = (9o,0,0, q), Pf, = {p ,psmd,0,pcosd), 
p'i-i = (p'o i P sm 4> sm Xj P sm 4> cos X, P cos </>) ■ 



C. Lorentz-invariant pair production and 
annihilation rates 

The equation for pair production and annihilation rate 
i? 77 4-> ee (_R12.i-i.34) in Eq. (|2ip is similar to the equation for 
the rate of muon production in photons fusion that we 
presented in [l|: 



R 



77~f4-ee 



2(2tt) 8 J 2E'{ 



d 3 pj f d 3 p1 f d 3 p% f d 3 p] 



222 



2ES 



2ES 



xS 4 ( P J 



Pi 



pI-p\)YMpI\ m i 



\pIpDY 



spin 



*fM)fM)fe{pl)h{p%)T- 2 T- 2 e^'l+^/ T , (35) 



Using Eqs. ([37)1 and delta functions we obtain the 
equations: 



pl-p 2 + -+ PoQo - pq cos (9) - m 2 = 0; (40) 

pl -P 2 + -7 -PoQo +pq cos {8) -to 2 = 0; (41) 

Po 2 ~P' 2 + 1+ Polo -p'gcos(^) = 0; (42) 

Po 2 -P ,2 + T -^0+^008(^=0; (43) 



then using S(f(x)) = J2i 1 /\f( x i)\$( x ~ x i)> we can 
rewrite the integral (similar to [6|) as 



8 



R 



77<->ee 



77. 7 61R / dq Q / dq I dpo / dp' / dp / dp/ / d(cos6») / d(cos0) / dxV|M 7 

^7TJ ID J 2 m JO Jqi J ql Jo Jo J -I J -1 JO 

(» - (>« + J - <nf°) * (✓ - (# + J) " 2 ) * (cos* - * 



gopo 



(44) 



where 



g / 4m 2 
,9 

?i,2 = ±2- 



(45) 
(46) 



qi,2 and gj 2 come from constrains cos#, cos</> < 1 and 6 
functions in Eq. (l4"4"l) . 

The integration over p, p 1 , cos 9, and cos </> can be done 
analytically considering 5 functions. Other integrals can 
be evaluated numerically. In order to simplify numerical 
integration we introduce dimensionless variables: 



q = (ql - Am 2 ) l ' 2 z- 



q / 4m 2 
p = -Sjl-—x; 



Pa = -jVi 



(47) 

(48) 
(49) 



< z < 1; -1 < x(y) < 1. 

In these variables, using the 8 function in Eg. (1441) . we 
obtain for u and t: 



u=(p-p) =m - - + 



Am 2 



(xy - \/(l- x 2 )(l - y 2 ) sin x) ; 



t = (p+p'Y = m z ----\jl- 



Am 2 



[xy - v/(l-:r 2 )(l-y 2 )sin X ). (50) 
Using new variables (|i7 |) -P5 ]l the Eq. J3S} is 



2 12 7T t > 



2ni 



dg exp(g ) / dzz 2 (q 2 - Am 2 ) 3/2 



f2lT nl *l / 

^ I d x I dx I dy\/l l M 



ii 



-l J-i 



77— ^ee I 



xT- 2 /7(|+P0)/7(|-P0 



: T-/ e (| + ,o)/e(f-po). 



(51) 



The matrix element is defined by Eq. [ 



D. Compton scattering Lorentz-invariant rate 



The equation for the Compton scattering rate is 



IK = l^/5f (52) 



(2tt) 8 J 2E% J 2EJ J 2E% 



xS 4 (pJ+Pt-pl-p^KpJpt \M ie ^ e \plpt) 



spin 



x/7(Pl7)/7(P27)/ e (P3)/4(P4)T- 2 T; 2 e -^+f7)/T 



The Feynman diagrams for Compton scattering are 
shown in figure Q] (b). The matrix element for Comp- 
ton scattering, Ea.([2"7]). can be obtained by crossing t 
and s, Eq. (l2~5l) . in the matrix element for pair production 
and annihilation, Eq. ([26|) . 



In this case s = (p3 + Pi) 2 , t = (pi — P2) 2 and u = 

(P3 -P2) 2 - 

We need to cross pi and P3 in variables g, p and p' 
definition, Eq. (|5o]l . 



V2 „f2 



In this case we have q = q$ — q = q' Q — q' 
and 



s > m 



Pi 



P2 



P3 



Pi 



(53) 



Then, similar to pair production and annihilation rate, 
we obtain: 
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R, 



■erf 



9eg~t 

(2tt) 8 



d\ / d*p / d±p'5 (p?) 5 (p 2 . - m 2 ) 6 (pi - m 2 ) <5 (p 2 ) (p°) (p°) 



(P§) 



x ]T |M 7e ^ 7e | 2 T- 2 / e (p°) / 7 (p°) T- 2 / 7 (p°) / e (p°) exp fo/T). 



(54) 



Using Eqs. ([55)1 and <5 functions from Eq. (|5"4"]) we ob- 
tain the equations: 



P' 2 -P' 2 



+ t - p' q +p'q cos {(f>) =0; (58) 



■ V + j + Po<7o - pq cos (0) = 0; 



(55) 



pl-p 2 + Paqa + pq cos (0) - to 2 = 0; (56) 

s then using S(f(x)) = J2i 1 /\f( x i)\ S ( x ~ x i)> we can 

p' Q 2 — p' 2 + — + PoQo — p'q cos (</>) — m 2 = 0; (57) rewrite the integral as 



2g e g 7 



m JO 



91 J H 
2\ V2\ 



oo poo 



rq2 rq 2 

Re 1 = ^^)6l 6 / dq ° I d<1 I dp ° I dp I d P I d P' I d ( C0S 9 ) I d ( COS & I d * 



JO 



^\M^ ei \ 2 5 \p-[p 2 + l-^) )S[ P '-[p^-^ + j) )S[oo S(P -^ + ^ 



2 4 



1/2N 



qoPo , m 



qp 1 2qp' 



x<5(cos0-B^-^-J/ e (|- p ) / 7 (| + p ) T c 7 2 T 7 - 2 / 7 (| - rt) /e (| +p>) exp fo/T), (59) 



where in this case 



91,2 



m 2 go ± Q( 1 m 
2s 2 



9l,2 



m 2 g q _ { 
2s 2 



to 
s 



(60) 
(61) 



The integration over p, p' , cos 9 and cos</> can be done 
analytically considering 6 functions. Other integrals can 
be evaluated numerically. In order to simplify numerical 
integration we introduce dimensionless variables: 



1 = (<7o 



m 2 )^ 2 z- 



Po 



m 2 q qf } _ m2 
2s 2 



, m 2 (?o }/, m 2 



(62) 
(63) 
(64) 



< z < 1; -1 < x{y) < 1. 

With these variables, using the (5 function in Eq. (l59"l) . 
we obtain for u and t: 



x(xy - V(i - x 2 )(l -y 2 )smx); 



n S TO S , 

= m 2 1 

2 2s 2 



f = (p + p 

x (a;y - v(l - x 2 )(l - y 2 ) sin x) 



rn 

s 



(65) 



Then limits for t are > t > — s + 2m 2 — m 4 / s. In the 
limit s >> m 2 the variables u and t for Compton scat- 
tering go to the same limit as u and t for pair production 
and annihilation, Eq. ([50|) . 

Using new variables ([62 ]) - (|6g]) the Eq. (J52j) is 
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9l9e 
2 n 7T 6 



dgo ex P(<7o) / dzz 2 (q 2 — m' 



xT: 



2 U (f "«,) A (I +Po) T-V, (f -/o) /e (f +P6 



(66) 



E. M0ller and Bhabha scatterings 



In our range of temperatures we assume: 



In this section we consider M0ller and Bhabha scat- 
terings, Eq.© and Eq.(J7J(, respectively. To remove 
Coulomb divergence in reaction rates we introduce ther- 
mal photon mass induced by Debye screening (l5j : 



(67) 



where lo p i is plasma frequency. 

Ae 2 f°° 

m 7 = T~2 / d PPfp(p)- 
Jo 

For highly relativistic plasma T >> m e 



y/Ar, 



-T. 



(69) 



The corresponding thermal mass of electron for high tem- 
peratures is jl 51 ] 



\J Ana 

m = r=-T. 

2V2 



(70) 



m = m„ + 



Ana 



(71) 



In this section m e = 0.5 MeV is electron mass without 
thermal effects. 



The Feynman diagrams for M0ller (a) and Bhabha 
scatterings (b) are shown in Fig. [2] and the matrix ele- 
ments are defined by Eas.(|29p and Eq. {J3DJ) , respectively. 
In both cases all initial and final particles are massive 
and have the same mass. After transformations similar 
to the ones used for Compton scattering the equation for 
the M0ller (Bhabha) rate is 



1 2g 2 f°° r^/il- s fte fi2 r°° r°° r 1 r 1 r 

Ree = —^L dq °L dq L dpo L dp ' o l dp l & }_«™°) I*™*) 1 



dx 



Y,\M ee \ 2 S[p- 



, 1/2 



(p'-(p'o 



- (Pa-™* + ~ 



1/2 



S cos 6 — 



QoPo 
qpf 



xS cos 9 — 



pq 



(72) 



9i,2 = Tt 



K d-—\ 



£ 4 /fi_^Y 



(73) 



(74) 



qi^2 and q\ 2 come from constraints cos 8, cos <f> < 1. Then 
we introduce dimensionless variables similar to Eqs. (ff)2")) - 



q = {q 2 a - Am 2 ) 1 ! 2 ?,- 

q If Am 2 \ 
P0= 2^{ 1 -—) X > 



Am 2 \ 



(75) 
(76) 

(77) 
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< z < 1; -1 < x(y) < 1. 

With these variables we obtain for u and i: 

« = {P-P) =2m 2 -- + - fl- — - 



x(xy- \J(l - x 2 )(l - y 2 )sinx); 

4m 2 



i = (p + p) = 2m 



s s 
2 ~ 2 



1 



x(xy— y/ (1 — x 2 )(l — y 2 ) sin%). 



In new variables ([75 |) - ([77|) the equation for scattering 
rate (El is 



(78) 



-Re 



1 2 



n.,6 



dq exp{q ) dzz 2 {q 2 - 4m 2 ) a/2 / <ix 

2m JO JO 



dtf 1- 



4m 2 



Po 



E 



|M ee | 2 / e (|-p )/ e (|+^)/ e (|-p 



(79) 



TABLE II: Values of N and A for the power law tempera- 
ture dependence of thermal cross section high energy tail, see 
Eq. (|34[) for the shape assumed, for cross sections presented in 
figure [6] 



Reaction 


N 


A [mb] 


7e <-> je 


-1.73 


120 


77 O e + e~ 


-1.73 


80 


e + e~ o 77 


-1.73 


46 



IV. NUMERICAL RESULTS AND DISCUSSION 
A. Thermal cross section and Thompson limit 

The thermally averaged cross sections (va) / c in ob- 
server frames for pair production, annihilation and 
Compton scattering Eq.(j4|), © are evaluated using the 
Lorentz- invariant rates, Eq.(|20|) with rates, evaluated nu- 
merically using Eq. ([5ip for pair production and annihi- 
lation and Eq. ([66|) for Compton scattering. In Fig. [6] 
we show (vcr) jc for Compton scattering (dashed red 
line), pair production (solid blue line) and pair anni- 
hilation (dash-dotted green line). For high T we ob- 
serve the power law fall-off of the thermal cross sections 
and the corresponding parameters are presented in Ta- 
ble [TTJ The power N — 1.73 in power law for pair pro- 
duction/annihilation and Compton scattering is close for 
nonthermal cross section in this case. At high tempera- 
tures these cross sections are almost parallel. This is in 
agreement with f-averaged cross sections Eqs. (l3"TT) - (|3"2")) . 

Considering electron production, the density of pho- 
tons with energy larger than the threshold for pair pro- 
duction drops in the tail of Boltzmann distribution with 
decreasing temperature, therefore the thermal cross sec- 



tion (va) starts to decrease at T < to and goes to as 
T — > 0. For electron-positron annihilation (va) jc stays 
for T — > 0, v — > 0, since a oc 1/v. 

At T << to, the Compton scattering cross section (va) 
is approaching the Thompson limit very slowly, which in 
case T = 1 is slightly [(va) comp (T = l)/ca th w 1.36] 
above the actual classical Thompson cross section 



87ra 2 
3to 2 



6.7 x 10 2 mb. 



This result is shown in Fig. [7] (dashed red line). This 
difference with the Thompson limit is due to quantum 
effects from massless photons. We can use Boltzmann 
distribution for photons only when T 7 << 1. In this case 
we obtain that thermal Compton cross section goes to the 
Thompson limit (green solid line in Fig. [7]) . In order to 
explain the difference between the cross section for T 7 = 
1 and the Thompson limit, we assume that electrons are 
at rest at small T and thermal photons are scattering 
without energy change JSj = Ef. Then the cross section 
quantum enhancement factor Q is the thermal average 
from the Bose enhancement factor / 7 (p) + 1: 



f-h 



27r 2 n 



f(p)(f(p) + l)p 2 dp. 



(80) 



7 JO 



1 and Q -¥ 1. For 
Further discussion 



In the Boltzmann limit f(p) + 1 - 
T 7 = 1 we obtained Q = 1.38 
about the Boltzmann and the quantum limit compari- 
son is found in Sec. IIVCI 



B. Mean Free Paths 

In Fig. |S] (left panel) we show the mean free paths 
of electron (positron) in EP 3 plasma in pair annihilation 
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1 T[MeV] 

FIG. 6: The thermally averaged cross sections (va) /c for 
Compton scattering (dashed red line), pair production (solid 
blue line) and annihilation ( dash-dotted green line) in the 
observer rest frame shown as a function of temperature T at 
T = 1 




FIG. 7: The thermally averaged cross sections (va) /c for 
Compton scattering for T 7 = 1 (dashed red line) and T 7 = 
0.01 (solid green line) as a function of temperature T. 



and Compton scattering reactions Eqs.Q and ([5]) (green 
dash-dotted and red dashed lines, respectively), and in 
M0llcr and Bhabha scattering reactions Eqs.© and 



mean free paths are evaluated according to Eq.([2T 



L 



ee— »77 



R, 



n e 



with rates evaluated, using Eq. (|51j) for pair annihilation, 
Eq. ffJfJ)) for Compton scattering, and Eg. ([79)1 for M0ller 
and Bhabha scatterings. To determine the possible ra- 
dius of opaque and equilibrated plasma [equivalent to 
opaqueness condition, Eq. ([2]), being satisfied] at a given 
energy, we show the radius of a plasma drop as a function 
of temperature [using Eq. (jll[) ] at expected experimental 
energies 500 J (upper light blue dotted line) and 50 J 
(lower purple dotted line). 

The mean free path of electron or positron in M0ller 
and Bhabha scattering remains much below the equilib- 
rium plasma radius at energy 500 J and the tempera- 
ture range considered here, up to the smallest temper- 
ature considered T = 0.3 MeV. Therefore electrons and 
positrons are thermally equilibrated at our temperature 
range. The L e in M0ller and Bhabha scattering is much 
smaller than L e in Compton scattering and in pair an- 
nihilation. Therefore electron and positron become ther- 
mally equilibrated first after their production by the laser 
field and then they produce photons over the larger time 
scale. 

We do not continue calculations for smaller tempera- 
tures since for (T << m e ) it becomes impossible to cre- 
ate equilibrium plasma with zero chemical potential con- 
sidered here, and hence for this reason our calculations 
of M0ller and Bhabha scattering rates are developed for 
temperatures T >> m. 

In Fig. [5] (right pannel) the mean free photon paths 
L 7 are shown for e + e~ pair production reaction (solid 
blue line) and Compton scattering (dashed red line) as 
functions of temperature T for plasma in chemical equi- 
librium (Tj = 1). Equation (l2~Tj) is used together with 
Eqs. (1511) and (|66[) for pair production and Compton 
scattering rates: 



R. 



L, 



R, 



The mean free paths for photon and electron scattering 
presented in Fig. [8] have power law distribution at high 
temperature, T >> m, 



L = A 



—) 

1 MeV J 



N 



(81) 



The coefficients A and ./V for mean free paths are listed 
in Table HIl 

We see that L 7e is less than L 77 at T = 1 over the 
full range of temperatures considered. The line for L ie 
follows that for L 77 at T > m e and L 7e is about 4 times 
smaller than L 77 , when values of L e7 are about twice as 
low as L 



at high temperatures. The factor making 
L ie to be half the length L ei is the photons scattering 



together (purple solid line) (T e 



T = 1). The on both electrons and positrons. 
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FIG. 8: Left panel: the mean free paths of electron (positron) shown as functions of plasma temperature in reactions Eqs.@, 
©(red dashed, green dash-dot lines, respectively), and reactions Eqs.© and © together (purple solid line); right panel: the 
mean free path of photon for Compton scattering and pair production at T = 1 (red dashed and blue solid lines) as functions 
of temperature T; radius of equilibrium (T = 1) plasma at energy 500 J (upper, light blue, dotted line) and at energy 50 J 
(lower, purple, dotted line) as a function of T on both sides. 



TABLE III: The values of N and A in the power law temper- 
ature dependence of free mean path high energy tail, Eq. ([8Tl 



reaction 


N 


A [nm] 


Lj, 7e 


- 1.3 


1.8 


L y , 77 •<-> e + e~ 


-1.3 


6.5 


L e , je o 7e 


- 1.3 


5 


L e , e + e~ o 77 


- 1.3 


2 


Lei 66 y 66 


- 1 


0.004 



TABLE IV: The values of maximum radius, minimum tem- 
perature and corresponding e + e~ pair and photon densities 
for opaque plasma at fixed energies 500 and 50 J 





opaque 


opaque 


opaque 


opaque 


E [J] 


Rmax [nm] 


T [MeV] 


n 1 [cm -3 ] 


n e [cm -3 ] 


500 


35 


0.5 


3.9 x 10 30 


5 x 10 30 


50 


3 


2 


2.55 x 10 32 


3.8 x 10 32 



At small temperatures, T < m, the photon mean free 
path L 7 drops with temperature faster than the plasma 
radius at a fixed energy. Plasma loses opaqueness when 
L 7e approaches the plasma radius in magnitude. For 
T m 0.5 MeV the free photon path L je is much smaller 
than radius of plasma, R = 35 nm, at energy 500 J. 
Since plasma is opaque, photons will stay in plasma vol- 
ume for a long time. For T < 0.5 MeV plasma becomes 
nonopaque and can lose energy, radiating photons. The 
plasma, which is opaque for photons, cannot have a ra- 
dius larger than R ss 35 nm for the considered energy 
content. In Table IIVI we show the values of maximum 
radius, minimum temperature and corresponding e + e~ 
pair and photon densities for opaque plasma at fixed en- 
ergies 500 and 50 J. 

Because at high temperatures the mean free paths de- 
crease slower than the plasma radius, R cx T 4 / 3 . There- 
fore there is a limit on the minimum energy need to cre- 
ate opaque plasma. A similar energy limit also exists for 



chemically equilibrated plasma, condition Eq.([3]). The 
energy E — 500 J is not large enough to satisfy this equi- 
libration condition exactly, but it is close to the minimum 
of required energy. The largest ratio i?/L 77 sa 2 arises 
at temperatures T w 4 — 6 MeV (Fig. [8j right panel) 
and plasma may be close to chemical equilibrium at this 
temperature range. The corresponding radius is ss 2 nm. 
The pair production or annihilation relaxation time Eq. 
(|16p is approximately 10 _3 fsatT = 5 MeV and the cor- 
responding plasma lifespan must be larger than that to 
satisfy the chemical equilibration condition. 

In this work we considered the highest density of chem- 
ical equilibrium plasma at T e = T 7 = T = 1 (for 
photon opaqueness condition T e = 1 is enough), be- 
cause the opaqueness condition Eq.® is quickly violated 
with Tj decreases. This is because the photon free path 
L 1 oc 1/T, when the radius of plasma R p \ oc 1/Y 1 / 3 , 
Eq. (fl~2|) . The result is i? p i/L 7 oc T 2 / 3 at a fixed energy 
content. Only for the relatively high plasma temperature 
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io- 1 10° 10 1 

T[MeV] 

FIG. 9: The photon mean free path for Compton scattering 
(lower red lines) and pair production (upper blue lines) ob- 
tained using reaction rate Eq. (|66[) (solid lines) and using rate 
from Eq. ([82]) (dashed lines), Boltzmann limit 

is it possible to have a below equilibrium density plasma 
with small T, without violation of opaqueness condition 
Eq.©. For example at T = 20 MeV and energy 500 J 
Rp\/Lry e f=a 8. In this case minimum possible T 0.3 
and radius at fixed plasma energy increases by factor of 
« 1.4 only (sa 0.4 nm), as the result R v \/L ie w 3. 

C. Lorentz-invariant rates and Boltzmann limit 

In order to evaluate the reaction rate in two - body pro- 
cesses in the relativistic Boltzmann (classical) limit, we 
can use reaction cross section a(s), and the relation [17] : 

Rim = i^TTT^ ds<r(s)^K,(VS/T), 

™ (82) 

(however, compared to Ref. [17| Eq.(17.16) iii2-^34 — > 
•ffi2->34/(TiT 2 )), where 

A 2 (s) = (s - (mi + m 2 f){s - (mi - m 2 ) 2 ), (83) 

mi and m 2 , gi and g 2 , Ti and T 2 are masses, degeneracy 
and fugacities of initial interacting particles. 

The cross sections can be evaluated using Eq. (|3Tj)- 
([33| . In cases considered here Eq. (|82|) with t-averaged 
cross section gives the same result for reaction rate as 
Eg. (122)) with T — > 0. In the general case the cross section 
averaging over t may also change the result compared to 
Eq. (f2"2"j) . Calculating mean free path L in the Boltzmann 
limit we also took the photon density in the Boltzmann 
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FIG. 10: The energies of plasma drop with different radius 
(R — 1/im, 100 nm, 10 nm, 1 nm and 0.1 nm) as a functions 
of T (dash-dotted lines - light blue, purple, brown, grey and 
blue, respectively); solid black line is energy vs. minimum T 
of opaque plasma; dashed blue line is energy vs. minimum T 
of chemically equilibrated plasma;) 

limit as n 7 — » n 7 /T at T — > to exclude all quantum 
distribution effects. 

In Fig. [S] we show our mean free paths for a photon 
subject to Compton or pair production process, for both 
quantum and Boltzmann (classical) limit as a functions 
of T. For Compton scattering quantum effects reduce the 
mean free path at high density (i.e. high T). Somewhat 
smaller, but a noticeable quantum effect remains at low 
temperature (T < m). For pair production there is no 
reduction of the mean free path at small temperature 
because of the threshold energy of photons in the pair 
production rate, and the smaller reduction for high T 
compared to Comptom scattering. 

V. SUMMARY AND CONCLUSION 

Current laser pulses can have energies up to 1 kJ and 
about 10 kJ is expected to be produced in the foresee- 
able future [Tf3]. Since we need to focus energy to sub 
diffraction size limit, R << A, only a part of the original 
energy is likely to be available to form plasma. 

The main purpose of this paper has been the evalua- 
tion of the domain of plasma size and temperature where 
thermally and chemically equilibrated EP 3 can exist at 
a given fireball energy. We evaluate the mean free path 
length L for photons and electrons in EP 3 plasma, assum- 
ing thermal equilibrium and using the Lorentz invariant 
reaction rate. This method allows us to take into account 
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quantum effects in a dense medium. 

Comparison of the electron (positron) mean free path 
with the plasma radius (figure El left) demonstrated that 
electrons and positrons can be thermally equilibrated by 
M0ller and Bhabha scattering at the temperature range 
that we considered. The small free path of electron and 
positron by M0ller and Bhabha scattering, compared to 
e + e~ pair annihilation, also means that electrons and 
positrons first become thermally equilibrated and then 
they produce thermal photons. 

Comparison of photon mean free path L 1 to plasma 
drop size R v \ (see figure [51 right) indicates that we may 
have restrictions on maximum plasma size, minimum 
temperature and minimum energy of opaque and chem- 
ically equilibrated plasma, defined by opaqueness con- 
dition Eq.©. We find that the mean free path for -fe 
Compton scattering is shorter than for pair production 
77 — > e + e~. Therefore opaqueness condition for photon 
will have L 7 = L ie ^ ie and differs from chemical equili- 
bration condition Eq.©. 

The practical summary of our results is presented in 
figure [TDJ The energies of different size plasma drops are 
shown as a function of temperature. The solid black line 
indicates the boundary below which the opaque plasma 
cannot exist; condition Eq.([2]) is violated. The blue 
dashed line indicates the boundary above which we have 
not only opaque but also chemically equilibrated plasma, 
condition Eq.Q. The green dotted line shows photon 
chemical freeze-out, L 1 fR v i = 1. The straight lines show 
constrains given by plasma size between E and T. 



In Fig. [TU] we see that at fixed plasma size to have 
opaque plasma less energy is needed than to have chemi- 
cally equilibrated plasma; for fixed energy the transition 
to chemically equilibrated plasma takes place at a smaller 
plasma sizes and higher temperatures and energies than 
the transition to opaque plasma. 

Only for E > 50 J can we expect to form an opaque 
plasma for T > 2 MeV. As more energy becomes avail- 
able, the temperature limit drops and the radius increases 
(R > 5 nm). It is interesting to note that much more en- 
ergy is needed to produce opaque plasma at lower com- 
pared to higher temperature, because plasma drop size 
has to increase considerably to satisfy the opaqueness 
condition Eq© . Thus one of main objectives of labora- 
tory experiments should be to seek paths to reduce the 
size of the domain in which energy deposition occurs. 

To summarize, our results indicate that chemically 
equilibrated, opaque, small - size is about 5-50 nm - 
plasma drop can be formed for energy of a few kJ which 
may be experimentally attainable in foreseeable future. 
Such a plasma drop would be an extraordinarily intense 
source of radiation and positrons. 
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